Quantum Zeno and anti-Zeno effects by indirect measurement with finite errors 
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We study the quantum Zeno effect and the anti-Zeno effect in the case of 'indirect' measurements, 
where a measuring apparatus does not act directly on an unstable system, for a realistic model with 
finite errors in the measurement. A general and simple formula for the decay rate of the unstable 
system under measurement is derived. In the case of a Lorentzian form factor, we calculate the 
full time evolutions of the decay rate, the response of the measuring apparatus, and the probability 
of errors in the measurement. It is shown that not only the response time but also the detection 
efficiency plays a crucial role. We present the prescription for observing the quantum Zeno and 
anti-Zeno effects, as well as the prescriptions for avoiding or calibrating these effects in general 
experiments. 
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It was predicted in a classic paper |Q that repeated 
measurements on a quantum unstable system, at time in- 
tervals r m , suppress the decay of the system for small r m 
— the so-called quantum Zeno effect (QZE) . It was as- 
sumed there that each measurement is completely ideal, 
i.e., it takes only an infinitesimal time, there is no error 
in the measurement, and the post-measurement state is 
exactly given by the projection postulate. However, any 
physical experiments do not satisfy all of these assump- 
tions, hence more careful studies have been desired. 

According to the general measurement theory, not only 
the unstable system in question but also a part of the 
measuring apparatus should be treated as a quantum sys- 
tem subject to the Schrodinger equation [|[ ||, f§ f§- It 
was clarified by such theories that the response time r r 
of the apparatus corresponds to r m of Ref. [Q. However, 
effects of the errors in the measurement are yet to be ex- 
plored, because the probability e of getting an erroneous 
result is determined not only by a finite response time r r , 
but also by the detection efficiency 1 — £oo (i.e., the ap- 
paratus occasionally fails to detect the decay even after 
an infinitely long waiting time). Moreover, these pioneer- 
ing theories, as well as pioneering experiments f?l 0, 
studied the case of 'direct' measurements, where the ap- 
paratus acts directly on the unstable system (e.g., shines 
laser light to excited atoms). In such a case, however, 
the dynamics of the unstable system would be affected by 
the apparatus even if the unstable system were a classical 
system, and thus the Zeno effect in direct measurements 
might not be peculiar to quantum systems. Hence, the 
most interesting case of 'indirect' measurements is yet to 
be explored, where the apparatus does not act directly 
on the unstable system, but detects a signal mediated by 
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some field. A promising theory of an indirect measure- 
ment was developed by Schulman j^] . However, since his 
model was an abstract one, application to real physical 
systems is not straightforward. For example, it did not 
clarify the conditions to observe acceleration of the decay 
by measurements, the anti-Zeno effect (AZE), which has 
also been attracting much attention B, [f(], [nj . To apply 
real experiments on QZE and AZE, more realistic models 
should be analyzed. Such analysis should also be impor- 
tant to general experiments, because the QZE or AZE 
might slip in advanced experiments in the near future. 
The purpose of this paper is to present a theory that sat- 
isfies all these requirements, using a realistic model of an 
indirect measurement with finite errors. 

The total system in our model is composed of three 
parts, (i) an unstable two-level system, which is initially 
in the excited state |x) with the transition energy f2o to 
the ground state |g), (ii) a field whose eigenmodes are la- 
beled by a wavevector k with any dimension, a quantum 
of which is emitted by the unstable system when it de- 
cays to |g), and (iii) a measuring apparatus that detects 
the emitted quantum by absorbing it, from which an ob- 
server gets to know the decay of the unstable system. A 
typical example for (i) is an excited atom, which emits 
a photon upon decay, and the photon is detected by a 
photodetector such as a photomultiplicr. Therefore, we 
hereafter call (i), (ii) and (iii) as an 'atom', 'photon' and 
'detector,' respectively, although the theory is applicable 
to other systems as well. In this model, neither a pro- 
jection operator nor the interaction Hamiltonian of the 
detector acts on the atom. Moreover, the measurement 
is of negative-result type, where no signal is detected un- 
til the atom decays: The QZE or AZE occurs just by 
waiting for the decay. 

The role of a detector is to convert a photon into 
other kinds of elementary excitations, which finally yield 
macroscopic signals after magnification processes, usu- 
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ally obeying classical mechanics. As a model of (the rel- 
evant part of) the detector, we assume elementary exci- 
tations (e.g, electron-hole pairs) with a continuous spec- 
trum, into which photons are converted. By taking the 
energy of |g) zero, the Hamiltonian of the system is taken 
as follows (with % = 1); 



H = ft |x)(x| +Hi +n 2 



Hi 



n 2 



dk 



(g k \x){g\b k + H.c.) + e k btb k 



dk duj 



k Ckuj + H.C. ) + UJC^Ckuj 



(1) 

(2) 
.(3) 



Here, b k is the annihilation operator for a photon with 
energy e k , g k represents the atom-photon coupling. The 
photonic dispersion relation can be nonlinear (e k 9^ |fe|) 
as in, say, photonic crystals. Every photon mode is lin- 
early coupled to a continuum of bosonic |T^ | elementary 
excitations, denoted by c kuJ , in the detector, with the 
coupling constant ^ kuJ . The commutation relations are 



-•>) C k', 



normalized as [b k , b' k ,] = 5(k — k') and [c kl 
S(k — k')S(u> — a/). The lifetime r k of a photon is de- 
termined by Cfcoj- We will show later that r k ~ r r , the 
response time of the detector. In most experiments, the 
detector does not cover the whole solid angle around the 
atom. When a photon is emitted in the uncovered direc- 
tion, it cannot be detected and has a long lifetime. As we 
will demonstrate later, we can encompass such realistic 
situations b y a llowing fc-dependence of r k . We here put 
(few = y/Vk @, which results in r k = (2nri k y 1 . 

The photon-detector part, Hi +H2, can be diagonal- 
ized in terms of the coupled- mode operator [ fl4| , which in 
this case is given by B k ^ = a k (p)b k + J duj /3 fc (/i, w)c fcttJ , 
where a k {p) = y/rjk/(fi - e k + m-qk) and f3 h ((J,,u) = 
Vk/ip- — £fe + 'wr}k)(n — ui + iS) + 6(fi - w). The com- 
mutation relations for Bkp, is given by [B kf i,B k , ,] = 
6(k — k r )5(oj — to'). Inversely, bk is expressed in terms of 
Bkfi as bk — J dfi a k ([i)B kfl . The Hamiltonian H can 
then be rewritten as follows, 



H = n |x)(x|+ / / 6kdn nBl M B klt 



dk d/i 



x)(g|B fcM + H.c. 



We further rewrite Ti using B^ that is defined by 

1 f ,, \/Vk9k 



dk- 



fJL-e k - nrq k 



-Bkn, 



where g^ satisfies 



l.9p| 2 = Jdk 



VkQk 



■(4) 



(5) 



(6) 



so that B^ is normalized as [B^, B* ,] = 5(/i — //). Then, 
the Hamiltonian finally reduces to 

n = r! |x)(x| + w 1 + 7Y 2 , (7) 

Hi = J dp [fe|x)<g|S„ + H.c.) + fiB^Bfj] , (8) 



where H2 consists of coupled-modes which do not interact 
with the atom. Now the physical meaning of B^ becomes 
apparent; it is part of the coupled- modes with energy fi 
that interact with the atom, whereas the other part is 
isolated. In terms of such operators, the Hamiltonian 
is expressed in the renormalized form, Eqs. (Q) and (^), 
where the atom is coupled to a single continuum of B^ 
with a coupling constant which is called the form fac- 
tor of interaction. The form factor under measurement is 
determined by Eq. @, from ek, gk, and r\k- It should be 
noted that there exists a sum rule J d/i|g M | 2 = J dfc|<?fc| 2 , 
which holds for any functional form of r/k ■ 

We first estimate the decay rate T under the measure- 
ment by a lowest-order perturbation in gk- We note that 
the straightforward application of Fermi's golden rule us- 
ing Tli as the interaction term gives a wrong result, be- 
cause strong effects of the detector are not involved. It 
is essential to use the renormalized form Tii as the in- 
teraction term. We then obtain a simple formula for the 
decay rate under the measurement; 



2tt / dk 



VkQk 



fio - £fc - i7T??fe 



(9) 



which should be compared with the free decay rate, Tq — 
2tt J dk\g k \ 2 5(fl — e k )- Note that Eq. @ includes non- 
perturbative effects of r\k . The formula clearly shows that 
the most important effect of the measurement (i.e., of fi- 
nite rjk) is to renormalize the form factor g^, and that the 
QZE (or AZE) occurs through the renormalization. Note 
that the formula is general, which holds for any forms of 
£fc, gk and r/k, and for any dimension of k. It is applica- 
ble, not only to spontaneous decay of an atom, but also to 
many other unstable systems if their Hamiltonian can be 
approximated by Eqs. ([!])- @ p2]. Moreover, the formula 
is also applicable to the case where the detector does not 
cover the full solid angle, yielding the detection efficiency 
< 1. In fact, suppose that only photons which are emit- 
ted in some solid angle Sd in the three-dimensional space 
are coupled to the detector, i.e., rj k = for (6>, </>) £ Sd, 
where k — (k sin 9 cos <j>, k sin 9 sin cf>, k cos 9). Then, Eq. 
(0) yields the simple formula; 



T = 2tt / k z dk 



dS 



ses d 



Vkgk 



fin 



+ 2tt / k 2 dk 



dS\g k \ 5(Cl -e h ) 



(10) 



S^Sd 



where dS = d cos 9d<p. 

Now we embody the above general results in an exam- 
ple in three dimension, in which (i) rjk — r\ = (2ttt)^ 1 
for (9,(j)) € Sd, whereas i]k = for (9,<j)) ^ Sd, (h) 
(k = k (we take c = 1), and (iii) gk takes the Lorentzian 
form after the angular integration; k 2 f Sf _ Rj dS \gk\ 2 ~ 



(1 - £oo ) 7 A 2 /[(fc - 

£oo7A7p - k a ) 2 



k ) 2 + A 2 ], and k 2 J s ^ s d d dS \g k \ 2 = 

f A 2 ], where ko ^S> A. In a case 
where gk depends only on k, is simply given by 
Soo = 1 — |Sd|/47r. In general cases, however, de- 
pends on both Sd and gk- It will turn out that 1 — 
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FIG. 1: The decay probability 1 — s(t) (solid curve), the prob- 
ability of getting an erroneous result e(t) (broken curve), and 
the probability of getting a detector response r(t) (dotted 
curve), when Qq — ko = 0, A = IOO7, r\ = 27, and e x — 0.2. 



corresponds to the detection efficiency, i.e., the probabil- 
ity of errors approaches as t — > 00. The meanings 
of 7 and A become transparent for Qq — ko, i.e., when 
the atomic transition energy coincides with the center 
of the Lorentzian. For A ^> 7, the decay rate of the 
atom is given by 2-717, while the transition from the ini- 
tial quadratic decrease to the exponential decrease in the 
survival probability occurs at t ~ 7] = 2/ A, which is 
called the 'jump time' H [Tl| . Using Eq. (|Io|), the renor- 
malized lowest-order decay rate is evaluated as 



2 7 T7e 00 A 2 



2tt 7 (1 -£qq)AA 

(fi - k ) 2 + A 2 (fi - fco) 2 + ■ 



r(r?,£oc), (11) 



where A = A + ixr\. This indicates that the effect of 
measurement on the decay dynamics become significant 
only for large rj satisfying ij > A, i.e., 7] > (27rn) _1 = 
t — T r , in accordance with the pervious studies g P . 

To see what is going on, we now calculate the 
temporal evolution of the wavefunction from the ini- 
tial state |x, 0,0). This can be pursued analyti- 
cally for the Lorentzian form factor. By putting 
\il>{t)) = e" m *|x, 0, 0) = /(t)|x, 0, 0)+ / dfc f h (t)\g, k, 0) + 
J J dkduj fkw(t)\gi 0, few), we calculate three probabili- 
ties; s(t) = \f(t)\ 2 (survival probability of the atom), 
e(i) = J dfc|/fc(/j)| 2 (probability that the atom has de- 
cayed but the emitted photon is not absorbed by the de- 
tector), and r(t) = J J dfedw|/fc w (/j)| 2 (probability that 
the emitted photon is absorbed). Assuming fast classi- 
cal magnification processes, we can interpret r(t) as the 
probability of getting a detector response, whereas e(t) 
is the probability that the detector reports an erroneous 
result. One of the advantages of the present theory is 
that all of these interesting quantities can be calculated. 
For example, s(t) is given by s(t) = 47r 2 |ci23e _lu ' - 
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FIG. 2: Plots of t" 1 lns(t). A = IOO7 and £00 = 0.2 in both 
figures. Qo — ko is taken (< A) in (a), and 200 7 (> A) 
in (b). The dotted lines show the renormalized lowest-order 
decay rate, Eq. (II lh . 



jA[(A +pnri)(u i - k Q ] 



C23ie 



-i(i! 2 t 



C3i2e IW 3'| 2 , where c ijk is given by c ijk 



(A — pirrf) AA] [(cji — Wf)(ui — 
and U j (j = 1,2,3) 
are the solutions of the cubic equation, (ujj — Oo)(^j ~ 
feo+iAK^-fco+iA) = -n-fAluj - k +i( P A+ (l~p)A)]. 
Figure □ plots 1 — s(t), e(t), r(t) for 77 = 27. At the ini- 
tial time stage (t <C t = (2irrj)^ 1 ), e(t) increases almost 
in parallel with 1 — s(t), and the detection probability 
r(t)[= 1 — s(t) — e(t)] remains almost zero. Around t ~ r, 
the emitted photon is gradually absorbed and r(t) starts 
to rise. Hence, the photon lifetime r can be regarded as 
the response time r r of the detector. The probabilities 
finally approaches the asymptotic values, 1 — s(t) — > 1, 
e(t) —> Soo, and r(t) — > 1 

To see the temporal behavior of s(t) more clearly, we 
have plotted i -1 In s(t) as a function of time in Fig. || for 
several different values of the parameters. At the begin- 
ning of the decay (t < 73), i _1 his(t) decreases linearly 
as < _1 lns(f) = —(J dk\gk\ 2 )t = —ir-fAt, for any value 
of fio — &o and for any values of the detector parame- 
ters r\ and £oo. Then, for t > tj, lns(i) approaches a 
constant value, which is well approximated by T(j], 
(dotted lines). These plots demonstrate that the decay 
dynamics is well described, except for the initial devi- 
ation, by the exponential decay with the renormalized 
lowest-order decay rate. In fact, we can show analyti- 
cally that formula ( |ll| ) is a good approximation to the 
asymptotic decay rate if 7 <C \ flo — k + iA| 2 /A. 

There is a remarkable difference between Figs. 0(a) and 
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on the dotted line, which is given by 



|n -k |/A 



FIG. 3: The phase diagram for the QZE and the AZE for a 
case of Lorentzian form factor (valid for any Soc ). The solid 
curve divides the QZE region and the AZE region. The dotted 
line shows the value of n at which the decay rate is maximized 
for each value of |f2o — ko\. 



(b), where different values of |r2 — fc |/A are employed. 
In case of |f2 — k \/ A = [see Fig. 0(a)], the decay rate 
decreases monotonously by increasing 77, i.e., the QZE 
occurs at any value of 7/. Contrarily, in case of |f2rj — 
&o|/A = 2 [Fig. |||(b)], the decay is enhanced for small 
f] (= 3O7), while it is suppressed for large 77 (= 2OO7). 
Thus, the AZE takes place for small 77. By analyzing 
Eq. ( pr] ) as a function of |Oo — &o| an d T], the 'phase 
diagram' discriminating the QZE and AZE is generated, 
which is shown in Fig. ||. The 'phase boundary' (solid 
curve) is given by 



Tj^ = p - ko f - A 2 ]/ttA, 



(12) 



on which the decay rate is not altered from the free rate 
To, while the decay rate takes the maximum value, 



r(7/ m >, goo ) 
r 



= £00 + (1 - £oc) 



Kin — kn 



2A\n -k \ 



- 1 



,(m) 



|n -*o| -A]/tt. 



(14) 



(13) 



We find that r]^ and n^ do not depend on e^. 

We finally discuss the significance of our results, for 
experiments on the QZE or AZE, and for general experi- 
ments. As mentioned earlier, indirect measurements are 
necessary, which have not been performed yet, for the 
complete experimental verification of the QZE and AZE. 
For such experiments, formula (^J) gives the necessary 
condition: T of Eq. (|9|) should significantly differ from 
the free decay rate Tq. In the Lorentzian case, this can be 
decomposed into the following conditions: (i) r r should 
be short enough; r r < Tj (which is a well-known condi- 
tion), (ii) T] should not be close to the phase boundary 
(|l2l), and (iii) should be so small that the first term 
of Eq. ([ll]) becomes dominant. Moreover, the QZE or 
AZE should be chosen according to the phase diagram, 
Fig. H| e.g., the AZE is most detectable on the dotted 
line. On the other hand, in general experiments, one 
usually wants to avoid the QZE and AZE in order to get 
correct results. Considering recent rapid progress of ex- 
perimental techniques and diversification of experimental 
objects, we expect that the QZE or AZE would slip in 
advanced experiments in the near future. To avoid the 
QZE and AZE, one must design the experimental setup 
to break at least one of the above conditions. For ex- 
ample, when performing an experiment with a high time 
resolution such that r r < Tj, then Eqs. ( |Io|) and jTl| ) 
suggest that £oo should be increased. If £00 cannot be 
increased to keep the sensitivity of such high-speed mea- 
surement, then one should adjust parameters in such a 
way that Eq. (|l2|) is satisfied, or, one should calibrate the 
observed value using our results, such as Eqs. ([H]) and 
(0), to obtain the free decay rate. 



[1] B. Misra and E. C. G. Sudarshan, J. Math. Phys. 18 756 
(1977). 

[2] K. Kraus, Found. Phys. 11 547 (1981). 
[3] E. Joos, Phys. Rev. D 29 1626 (1984). 
[4] A. Beige and G. C. Hegerfeldt, Phys. Rev. A 53 53 
(1996). 

[5] E. Mihokova, S. Pascazio, L. S. Schulman, Phys. Rev. A 
56 25 (1997). 

[6] W. M. Itano, D. J. Heinzen, J. J. Bollinger, and D. J. 

Wineland, Phys. Rev. A 41 2295 (1990). 
[7] Chr. Balzer, R. Huesmann, W. Neuhauser and P. E. 

Toschek, Opt. Commun. 180 115 (2000). 
[8] M. C. Fischer, B. Gutierrez-Medina, and M. G. Raizen, 



Phys. Rev. Lett. 87 040402 (2001). 
[9] L. S. Schulman, Phys. Rev. A 57 1509 (1998). 
[10] A. G. Kofman and G. Kurizki, Nature 405 546 (2000). 
[11] P. Facchi, H. Nakazato, and S. Pascazio, Phys. Rev. Lett. 

86 2699 (2001). 
[12] In general, excitation of fermion pairs at low density are 

well described by a collection of harmonic oscillators. 
[13] It is straightforward to include u dependence of . We 

find that it just introduces uninteresting complexity. 
[14] K. Koshino and A. Shimizu, Phys. Rev. A 53 4468 

(1996). 



